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ABSTRACT 

Aims. Several frequency combinations are widely used in the analysis of stellar oscillations for comparisons between models and observations. 
In particular, the "small separations" can help constraining the stellar cores. We showed in a previous paper that they can change sign, in 
contradiction with the "asymptotic theory", and that this behaviour could correspond to signatures of convective and/or helium cores. Here we 
analyse this behaviour in detail by systematic modelling during stellar evolution. 

Methods. We computed evolutionary tracks for models with various masses (from 1.05 to 1.25 M Q ) and various chemical compositions, 
with and without overshooting. We computed the adiabatic oscillation frequencies of the models and analysed the evolution of their small 
separations along an evolutionary track. 

Results. We find that, for all cases, the stars go through a stage during their evolution, where the small separations computed between degrees 
£ = and 1 = 2 become negative in the observed range of frequencies. This behaviour is clearly related to the signature of a helium-rich core. 
We discuss the consequences for interpreting of the acoustic frequencies observed in solar-type stars. 



1. Introduction 

Asteroseismology has proved to be an excellent tool for deriv- 
ing stellar parameters more precisely than ever before: mass, 
age, radius, metallicity, helium abundance value, stellar grav- 
ity, effective temperature. This has been discussed recently in 
detail for the case of the exoplanet-host star t Hor (Vauclair et 
al. [2008 ). Spectroscopic observations lead to values of T e ff, 
g, and [Fe/H], which differ for different authors, according to 
the scale they use for the effective temperature determinations. 
The associated error bars are quite large. In this framework, ob- 
serving the acoustic oscillations of the star and determining the 
mode frequencies is extremely helpful. 

Once the oscillations of the star have been observed long 
enough (typically eight nights with HARPS), the analysis of the 
Fourier transform may lead to the determination of the "large 
separations", defined as Av,,/ = v„+i/ - v„/ and to the iden- 
tification of the modes, as discussed in Bouchy et al. (2005). 
Meanwhile, various evolutionary tracks are computed with dif- 
ferent input parameters (mass, chemical composition, presence 
or not of overshooting, etc). For a given set of parameters, 
only one model may reproduce the observed frequencies sat- 
isfactorily (Bazot et al. 2005; Soriano et al. l2007l Laymand & 
Vauclair 2007, Vauclair et al. I2008l l. As shown in Vauclair et 
al. ((2008), the various models obtained in this way have sim- 
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ilar gravities and ages. The other parameters are constrained 
with the help of the spectroscopic observational boxes. 

Asteroseismology can also give information about the in- 
ternal structure of the stars and, more specifically, about the 
regions where the sound velocity changes rapidly. This hap- 
pens in various transition layers, such as the limits of convec- 
tive regions or layers with strong helium gradients. The transi- 
tion layers that occur in stellar outer regions (bottom of outer 
convective zones, diffusion-induced helium gradients) may be 
characterised owing to the reflexion of the acoustic waves on 
the stellar surface and on the region of rapidly varying sound 
velocity. This is generally studied with the help of the "sec- 
ond differences", defined as S^v^t = v n +i,e + v n -\,t ~ 2v„/, 
which present oscillations with a period related to the acous- 
tic depth of the transition layer (e.g. Gough |T990l Monteiro & 
Thompson [T998l Vauclair & Theado I20041 Castro & Vauclair 
2006). On the other hand, the transition layers in the stellar 
internal regions are better characterised using the "small sep- 
arations", defined as dv n j = v n ,c - Vn-\,c+2 (e-g- Roxburgh & 
Vorontsov |1994| l. In a series of papers, Roxburgh & Vorontsov 
(see Roxburgh & Vorontsov 2007 and references therein) show 
how the presence of a convective core could lead to oscillations 
in the small separations. 

In a previous paper dedicated to the study of the exoplanet- 
host star HD 52265 (Soriano et al. 2007), we showed that in 
some cases the small separations, which should be positive in 
first approximation, could become negative. We explained how 
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this special behaviour was related to either a convective core or 
a helium core with abrupt frontiers, resulting from the presence 
of a convective core in the past history of the star. 

In the present paper, we give a general analysis of this effect 
for solar-type stars, with or without overshooting. We focused 
our analysis on stars that have either a solar metallicity or that 
are overmetallic, which is the case of exoplanet-host stars. We 
show that, along their evolutionary tracks, all stars go through 
a stage where the small separations between the £ — and 
I — 2 modes become negative, whether near the end of the 
main-sequence or at the beginning of the subgiant branch. This 
leads to a special behaviour in the "echelle diagram", which 
is obtained by plotting the mode frequencies in ordinates and 
the same frequencies modulo the large separations in abscis- 
sae (e.g. Soriano et al. 2007 ). Instead of remaining parallel, as 
predicted by the asymptotic theory (Tassoul 1980), a crossing 
point appears between the t — and 1-2 curves. This may 
also happen between the i = 1 and { = 3 curves, but it is less 
frequent and we concentrate here on the I — and t = 2 case. 

For comparisons with observations, this behaviour be- 
comes important if the crossing points correspond to frequen- 
cies below the cut-off frequency of the acoustic waves, which is 
the highest frequency for which the waves are reflected on the 
atmospheric regions. As this cut-off frequency depends on the 
characteristics of the stellar outer regions and may vary from 
star to star, we chose to study all the cases for which the cross- 
ing point occurs below 3.5 mHz. Along the evolutionary track, 
we computed the models for which this behaviour appears for 
the first time, for various masses, and give their ages. In the 
following, we refer to these special models as "transition mod- 
els". As the apparition of this behaviour is related to the radius 
of the convective or helium core, negative small separations ap- 
pear earlier when overshooting is introduced. 

This particular behaviour of the small separations could be 
used to derive the size of stellar cores and the presence of over- 
shooting for stars at the end of the main-sequence. A direct 
application to the exoplanet-host star /j Arae will be given in a 
forthcoming paper. 

2. Evolutionary tracks and models 

We computed series of evolutionary tracks using the Toulouse 
Geneva Evolution Code (see Hui Bon Hoa 2007 for a general 
description of this code) with the OPAL equation of state and 
opacities (Rogers & Nayfonov 2002; Iglesias & Rogers 1996) 
and the NACRE nuclear reaction rates (Angulo et al. 1999)). In 
all our models, we included microscopic diffusion as described 
in Michaud et al. (120041 ) and Paquette et al. ( 119861 ). 

The convection was treated in the framework of the mixing 
length theory and the mixing length parameter was adjusted as 
in the Sun: a — 1.8 (Richard et al. 120041 ). Models were com- 
puted for a range of masses from 1.05 to 1.25 M Q . 

We computed six series of models. The first three series 
were computed without overshooting and the second three se- 
ries with overshooting at the limit of the stellar core. Here over- 
shooting is described as an extension of the central convective 
zone by a length a ov Hp, where Hp is the pressure height scale, 
and a ov the overshooting parameter. In our computations, we 



fixed a ov to 0.20. In each case, the three series differ from their 
abundances: 

- Metallic and helium solar values. For the metallic values, 
the "old" Grevesse & Noels (1993 ) abundances, compati- 
ble with helioseismology, are used: X ml = 0.7097, T,„, 
0.2714. Z ini = 0.0189 (Fig. Q] upper panel, and Table Q] 
for the case without overshooting; Fig. [2] upper panel, and 
Table|4]for the case with overshooting). 

- Overmetallic abundance [Fe/H]=0.30, and helium value 
computed as for the chemical evolution of galaxies, namely 
dY/dZ = 2.8 ± 0.5: X M = 0.6 648, Y ini = 0.3027, Z M = 
0.0324 (Izotov & Thuan 12004] ) (Fig.[T] middle panel, and 
Table |2] for the case without overshooting; Fig. [2] middle 
panel, and Table|5]for the case with overshooting). 

- Overmetallic abundance [Fe/H]=0.30, and the solar he- 
lium value Y : X M = 0.6961, Y ini = 0.2714, Z M = 
0.0325 (Fig.Q] lower panel, and Table|3]for the case with- 
out overshooting; Fig. [2] lower panel, and Table [6] for the 
case with overshooting). 

We computed adiabatic oscillation frequencies for a large 
number of models along each evolutionary track, using the 
PULSE code (Brassard et al. 1 1992b . The oscillation frequen- 
cies were computed for degrees I = to I = 3, which are 
the only degrees observable for solar-like stars, because of the 
lack of spatial resolution. We only kept frequencies between 
1.5 and 3.5 mHz, corresponding to the typical observational 
range for solar- type stars. Their radial orders range between 4 
to 100. Using these frequencies, we computed the small sepa- 
rations for each model and analysed in what conditions these 
quantities can become negative. For all evolutionary tracks, we 
found models where the small separations become negative at 
some frequencies, which means that the I — and I — 2 lines 
cross over at a given point in the echelle diagram (cf Fig. [3] for 
example). In all cases, the frequencies of these crossing points 
decrease for increasing stellar age. For each computed track, 
we picked up the "transition model" for which the frequency 
of the crossing point is 3.5 mHz. These models are represented 
by crosses in Figs. [T] and [2 and their characteristics are given 
in TablesQ]to[6] 

3. Analysis of the results 

A first analysis of the possibilities of negative small separa- 
tions in stars was given in Soriano et al. (12007 ) for the specific 
case of the exoplanet-host star HD 52265. When we computed 
models for this star, taking spectroscopic constraints into ac- 
count, we were surprised to find two models that showed this 
specific behaviour: one model at the end of the main sequence 
with a mass of 1.31 M , which had a convective core, and one 
model at the beginning of the subgiant branch with a mass of 
1 .20 M , in which the convective core, present during the main 
sequence, had disappeared. For both cases we realised that the 
main reason for the negative small separations was related to 
the high helium content of the cores. The basic role of convec- 
tion in this respect was its action of concentrating helium inside 
a sharp core during the main sequence. 



Soriano & Vauclair: Signatures of convective and/or helium cores. 



3 




3.80 3.78 3.76 
log(Teff) 



cn 4.2 
o 

4.4 



. [Fe/H] = d.3b, V ' 








S - s * f 








L 1 ' ; ' f 

1.25 v > '■ i 




: i - 2 ° s --i.-i5".. V 


No \ 


1J0 1.05 


ov . 



3.78 3.76 3.74 
log(Teff) 





3.8 


. |:e/--| 0.30, Y 




4.0 




cn 






CP 

o 


4.2 


■ **** ^ X 




4.4 


1.25 v -,2 b -\ _ ' ■ no : 

U5 i.io \__ OV . 

1-05 






3.78 3.76 3.74 



log(Teff) 



Fig. 1. Evolutionary tracks computed for solar metallicity [Fe/H] (upper panel), overmetallicity [Fe/H]=0.30 with a helium 
abundance following the law of the chemical evolution of galaxies Y c (middle panel), and overmetallicity with a solar helium 
abundance Y Q (lower panel). All these tracks are computed without overshooting. The masses represented are: 1.05, 1.10, 1.15, 
1.20, and 1.25 M . The crosses represent the "transition model" for each evolutionary track, or first model from which we found 
negative small separations below 3.5 mHz. The characteristics of these models are displayed in TablesQ~|[2l and[3] 
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Fig. 2. Evolutionary tracks computed with the same hypothesis as in Fig.Q]but with overshooting. Here, the overshooting pa- 
rameter is fixed as a ov = 0.20. The masses represented are the same ones as in Fig. [TJ On each evolutionary track, we have 
represented the "transition model" by crosses, or first model from which we have negative small separations below 3.5 mHz. The 
characteristics of these models are displayed in Tables HJEI and[6] 
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Table 1. Parameters of the "transition models" computed without overshooting, with solar metallicity [Fe/H] Q and solar helium 
value (Fig.Q] upper panel). 



M*(Ms) 


Age (Gyr) 


x c 


Y c 


R*(cm) 


U/L 


log T c// 




R„/R 


tec (S) 


1.05 


7.3 


< 10~ 4 


0.9758 


9.07el0 


1.854 


3.771 


4.332 


(0.07) <n) 


(124.4) <fc) 


1.10 


6.0 


< 10~ 4 


0.9764 


9.54el0 


1.913 


3.779 


4.268 


(0.07) (a) 


(128.6) (i,) 


1.15 


4.9 


< 10~ 4 


0.9785 


10.12el0 


2.681 


3.788 


4.176 


(0.06) <a) 


(127)® 


1.20 


3.7 


0.1462 


0.8338 


9.87el0 


2.579 


3.789 


4.216 


0.05 


113 


1.25 


3.0 


0.1986 


0.7817 


10.30el0 


2.992 


3.796 


4.196 


0.06 


123.7 



<a) The values in parenthesis correspond to the cases where the convective core has disappeared, leaving a helium core with abrupt limits; here 

the values correspond to the linear radii of the helium cores. 

(i) The values in parenthesis correspond to the acoustic radii of the helium cores. 

Table 2. Same as TableQ]for the overmetallic models ([Fe/H]=0.30) with a helium mass fraction Yc (Fig.Q] middle panel). 



M*(Mo) 


Age (Gyr) 




Y c 


R.(cm) 


UIU> 


log T e// 


l°gg 




fee (s) 


1.05 


7.6 


< 10- 4 


0.9567 


8.75el0 


1.511 


3.757 


4.263 


(0.07) 


(121.6) 


1.10 


5.7 


0.0462 


0.9173 


8.79el0 


1.594 


3.761 


4.279 


0.05 


98.2 


1.15 


4.7 


0.1058 


0.8588 


9.15el0 


1.819 


3.767 


4.263 


0.06 


113.1 


1.20 


3.7 


0.1760 


0.7893 


9.50el0 


2.117 


3.776 


4.250 


0.06 


126.4 


1.25 


3.1 


0.2204 


0.7451 


9.98el0 


2.481 


3.782 


4.224 


0.07 


140 



Table 3. Same as TableQ]for the overmetallic models ([Fe/H]=0.30) with a helium mass fraction Y (Fig.Q] lower panel) 



M,(M Q ) 


Age (Gyr) 


X c 


Y c 


R*(cm) 


L*/L© 


log T c// 


l°gg 


Rir/R 


tec (S) 


1.05 


9.7 


< 10~ 4 


0.9545 


8.52el0 


1.294 


3.746 


4.285 


(0.06) 


(104.6) 


1.10 


7.9 


< 10~ 4 


0.9565 


8.94el0 


1.536 


3.754 


4.264 


(0.06) 


(108.1) 


1.15 


6.6 


< io~ 4 


0.9609 


9.52el0 


1.877 


3.762 


4.229 


(0.06) 


(111.4) 


1.20 


4.9 


0.1220 


0.8427 


9.30el0 


1.848 


3.765 


4.267 


0.05 


108 


1.25 


4.0 


0.1787 


0.7864 


9.70el0 


2.146 


3.773 


4.249 


0.06 


122.5 



We showed how negative small separations could be a signa- 
ture of the size of a helium core, as well as that of a convective 
core when it is still present. 

We give below a more complete analysis of the results we 
obtained by doing systematic studies of this effect for solar- 
type stars. We first recall useful theoretical points, then we dis- 
cuss our computational results. 

3.1. Theoretical analysis 

In all our models, the acoustic frequencies were computed pre- 
cisely using the PULSE adiabatic code (Brassard et al. 1992). 
The results take the particular features of the stellar interiors 
into account. That the small separations become negative in 
some cases is real. 

This seems surprising at first sight because it contradicts 
what is usually called "asymptotic theory", as developed by 
Tassoul (1980). According to this analytical description of the 
oscillations, the large separations should be constant, equal 
to half of the inverse of the acoustic time, defined as the 
time needed for the acoustic waves to cross the stellar radius. 
Meanwhile, the small separation should vary quite slowly and 



remain positive. Although these approximate expressions are 
not used in real computations, they are very useful for under- 
standing the underlying physics. 

As discussed in Soriano et al. (120071 1. the asymptotic the- 
ory has been derived with some assumptions that are no longer 
valid in the cases we are studying. The most drastic of these 
assumptions is related to integrals that should be computed be- 
tween the internal turning point of the waves r, and the stellar 
surface, whereas they are replaced by integrations from zero 
to R. This means that the sound velocity at the turning point 
is assumed identical to that in the stellar centre. This assump- 
tion is completely wrong in stars with convective cores, and 
it becomes worse and worse as the helium-to-hydrogen ratio 
increases. A second important assumption is that the large sep- 
arations for all frequencies and all modes are assumed identical 
to 



Avn = 



Jo c 



(1) 



This is not correct, as the ( + modes do not travel down to the 
stellar centre. They are trapped at the turning point r,, so that 
for each mode, integral (Q]) should be computed from r, to R, 
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Table 4. Parameters of the "transition models" computed with overshooting, with solar metallicity [Fe/H] Q and solar helium 
value (Fig. [2] upper panel) 



M*(Ms) 


Age (Gyr) 


x c 


Y c 


R*(cm) 


U/L 


log T e// 






tec (S) 


1.05 


7.3 


< 10~ 4 


0.9758 


9.07el0 


1.854 


3.771 


4.332 


(0.7) (o) 


(124.4) (6) 


1.10 


5.5 


0.2285 


0.7514 


9.02el0 


1.896 


3.775 


4.256 


0.06 


116.8 


1.15 


4.8 


0.2369 


0.7430 


9.67el0 


2.284 


3.780 


4.215 


0.07 


129.7 


1.20 


4.2 


0.2429 


0.7371 


10.28el0 


2.723 


3.784 


4.174 


0.07 


139.9 


1.25 


3.7 


0.2472 


0.7329 


11.04el0 


3.223 


3.789 


4.136 


0.07 


148.1 



<a) The values in parenthesis correspond to the cases where the convective core has disappeared, leaving a helium core with abrupt limits; here 

the values correspond to the linear radii of the helium cores. 

(i) The values in parenthesis correspond to the acoustic radii of the helium cores. 

Table 5. Same as Table[4]fbr the overmetallic models ([Fe/H]=0.30) with a helium mass fraction Yg (FigfJ] middle panel). 



M*(Mo) 


Age (Gyr) 




Y c 


R*(cm) 


U/L 


log T eff 






tec (S) 


1.05 


7.6 


< 10- 4 


0.9567 


8.77el0 


1.516 


3.757 


4.261 


(0.07) 


(121.6) 


1.10 


6.3 


0.2107 


0.7539 


9.05el0 


1.594 


3.757 


4.258 


0.07 


136.4 


1.15 


5.4 


0.2179 


0.7468 


9.58el0 


1.912 


3.763 


4.223 


0.07 


146.5 


1.20 


4.7 


0.2232 


0.7418 


10.20el0 


2.277 


3.769 


4.189 


0.07 


155.4 


1.25 


4.1 


0.2293 


0.7358 


10.83el0 


2.691 


3.773 


4.153 


0.07 


161.3 



Table 6. Same as Table|4]for the overmetallic models ([Fe/H]=0.30) with a helium mass fraction Y (Fig.[2] lower panel). 



M,(M ) 


Age (Gyr) 


X c 


Y c 


R*(cm) 


U/L 


log T e// 


logs 


R„./R 


tec (S) 


1.05 


9.7 


< 10~ 4 


0.9545 


8.54el0 


1.298 


3.746 


4.284 


(0.06) 


(104.4) 


1.10 


8.0 


< 10~ 4 


0.9565 


8.97el0 


1.543 


3.753 


4.261 


(0.06) 


(107.7) 


1.15 


6.6 


0.2267 


0.7379 


9.20el0 


1.627 


3.754 


4.259 


0.07 


136.3 


1.20 


5.7 


0.2347 


0.7301 


9.75el0 


1.938 


3.761 


4.227 


0.07 


144.8 


1.25 


5.0 


0.2402 


0.7248 


10.35elO 


2.290 


3.766 


4.192 


0.07 


155.6 



not from to R. As a consequence, the mean large separation 
for a mode £ is larger than Avq, due to the smaller integral. This 
effect can be important if the sound velocity decreases strongly 
in the central regions. 

Following Tassoul ( 1980), but relaxing these assumptions, 
Soriano et al. (2007) derived the following approximate expres- 
sions for the ^=0-^=2 and for the £ — 1 - i — 3 small 
separations: 



1 6Av 
6v 02 -(« + -+ a)(Av - Av 2 ) - — -= — '■ 

4 47T / V„_ 



Av 2 r R l dc 

Vn-1,2 J r , r dr 
Av 3 \ f R 1 dc 
Vn-1,3 / Jr, r dr 



3 \/ a an/ Ay l 6Av 3 

8v 13 * (n + - + ar)(Av! - Av 3 ) + — — 

4 \ln L v n ^\ 2n A v n - 

where n is the radial order of the modes, a corresponds to a 
surface phase shift, and Avo, Avi, Av2, and Av 3 are the large 
separations computed respectively for the degrees ( — 0, 1,2, 
and 3. 

From now on, we basically concentrate on the <5vo2 separa- 
tions, which are the most relevant for our purpose. We can see 
that the integral I(r) = J r \ j- r dr plays an important role in the 



computations of the small separations. Using the real boundary 
value r, instead of zero can lead to significant changes in the 
results. Also the difference between the large separations for 
the £ = and the I — 2 modes, which vanishes in the asymtotic 
theory, can be important, even if it is small. Indeed this differ- 
ence is multiplied by the radial order n which can be significant 
for the considered modes (of order 30). 

We will see below that for some models, the integral I(r), 
which is generally negative, may become positive at the bound- 
ary of the core when the helium over hydrogen ratio is high 
t^f^ygh. This is due to the resulting rapid variation of the sound 
velocity, which is significantly lower in the core than expected 
from simple extrapolation below the outer layers. In this case, 
*tfre ^ign inversion in the integral is sufficient to explain why the 
small separations become negative (Fig.|4|i. The frequency for 
which this occurs corresponds to that of the 1-2 waves for 
which the turning point is at the core boundary. 

In other cases, it may happen that the integral value changes 
abruptly at the core boundary, but not enough to become pos- 
itive. Then, a slight difference between the large separations 
for the £ — and the £ = 2 modes may be enough to lead to 
the inversion of the sign of the small separations (Fig. [3). In 
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frequency (mHz) frequency modulo 89/zHz 




Fig. 3. Small separations (upper-left panel), echelle diagram (upper-right panel), large separations (lower-left panel), and integral 
I{r) (lower-right panel) for a model with a solar metallicity, 1 .20 M Q and 3.683 Gyr. It is a main sequence model with a convective 
core. Here the integral I(r), although rapidly changing near the core, does not change sign. On the other hand, as can be seen in 
the lower-left graph, the large separations are higher for 1-2 than for t — 0, which is enough to induce a sign change in the small 
separations. For the echelle diagram, solid lines are for i — 0, dotted lines for 6=1, dashed lines for t = 2, and dotted-dashed 
lines for I = 3. 



this case, the frequency where this behaviour occurs may be 
slightly different from that related to the core boundary. 

As a conclusion, that the small separations change sign is 
related to the core boundary and may be physically understood. 
However the use of this sign inversion to derive the core size 
can only be done with precise model computations. 

3.2. Analysis of the small separations along an 
evolutionary track 

Figures [5]to[7]present the internal parameters for three models 
along an evolutionary track of 1.15 M with [Fe/H]=0.30, with- 
out overshooting. The first model lies on the main sequence, in 
a situation where the small separations are always positive in 
the considered frequency range (1.782 Gyr, Fig. |5). The sec- 
ond model is more evolved, and the small separations for the 
degrees I = - I = 2 become negative at a frequency of 3.4 
mHz (4.654 Gyr, Fig|6]l. It corresponds to what we have called 
the "transition model". The third model is still more evolved on 
the subgiant branch (5.937 Gyr, Fig. [7]). 

In each figure, four different graphs are presented: the small 
separations (upper-left panel), the echelle diagram (upper-right 



panel), the sound speed profile (lower-left panel), and the inte- 
gral I(r) (lower-right panel). 

At the beginning of the main sequence (model presented in 
Fig|5]l, the convective core of the star is still small (R cc /R = 
0.03). The central helium abundance is low and there is no dis- 
continuity in the sound speed profile at the centre of the star. 
The integral I{r) is negative in the whole star. In that case, the 
small separations <5vo2 are positive for the full range of frequen- 
cies considered, and there is no crossing point in the echelle 
diagram. 

When the age of the star increases, the values of the small 
separations decrease. The convective core develops, and there 
is more helium in the centre of the star. A discontinuity in 
the chemical composition and in the sound speed profile ap- 
pears. The model shown in Fig. [6] is the first model for which 
we found negative small separations below 3.5 mHz. We call 
it the "transition model" and show its position as a cross on 
the corresponding track in Fig.Q] The convective core is larger 
iXcclR = 0.06) than for the previous model. We can see the dis- 
continuity in the sound speed profile (Fig. [6] upper-left panel). 
As a consequence, the integral I{r) changes drastically between 
r = and r = r t = R cc . Here it does not change its sign but, 
as discussed before, its variations are large enough to lead to 
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frequency (mHz) frequency modulo 78/zHz 




Fig. 4. Small separations (upper-left panel), echelle diagram (upper-right panel), large separations (lower-left panel), and integral 
I{r) (lower-right panel) for a model with a solar metallicity, 1 .20 M and 4.280 Gyr. It is an evolved model (subgiant branch) with 
a helium core. Here the integral I(r) changes its sign near the core. This explains the sign change in the small separations that 
occur when the 1-2 waves reach the core boundary and the crossing point in the echelle diagram. For the echelle diagram, solid 
lines are for i — 0, dotted lines for t = 1, dashed lines for 6 = 2, and dotted-dashed lines for I = 3. 



negative small separations. The lines £ = and £ = 2 cross in 
the echelle diagram. 

In the third model presented in Fig. [7] the convective core 
has disappeared, leaving a helium core. There is still a sharp 
discontinuity in the sound speed profile at the edge of the core 
(Fig-El lower-left panel), and the small separations change sign 
at a lower frequency than for the previous model (Fig. [71 upper- 
left panel). 

The frequencies at which the small separations become 
negative go on decreasing while the star evolves along the sub- 
giant branch, until the stellar structure begins to change deeply, 
with an increase in the outer convective layers and a shrinking 
of the core. A discussion of this evolutionary stage is beyond 
the scope of the present paper and will be given elsewhere. 
Here we concentrate on the main sequence and early subgiant 
phases. 

3.3. Influence of the stellar mass 

Tables [T]to|6]present the parameters for the "transition models" 
for various physical assumptions. Let us recall the meaning of 
these models. As discussed in the previous section, the frequen- 
cies for which the small separations 6vq2 change sign decrease 



with increasing age along an evolutionary track. The transition 
model is arbitrarily defined on each track as the one for which 
this characteristic frequency is about 3.5 mHz. The interest of 
comparing such specific models is that it helps in analysing the 
influence of parameters like stellar mass or chemical composi- 
tion on the behaviour of the small separations. 

In each table, the results are presented for the same five 
masses, between 1.05 and 1.25 M Q . We can see that the ages 
of the transition models decrease for increasing masses. This is 
not only due to the more rapid evolution time scales, as the cor- 
responding helium mass fractions also decrease with increasing 
mass. On the other hand, we can see that the total radii of the 
models, as well as the radii of their convective or helium cores, 
do not show strong variations with mass. 

These results are quite understandable if we remember the 
reason for the appearance of the negative small separations. As 
discussed previously, they generally appear for frequencies for 
which the turning point of the t = 2 waves corresponds to the 
edge of the core. The relation between the frequencies of the 
waves and their turning points is given by (see Soriano et al. 
120071 

2nv = — yf£{fi + 1), (4) 
r 
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Fig. 5. Small separations (upper-left panel), echelle diagram (upper-right panel), sound speed profile (lower-left panel), and 
integral I(r) = £ \j- r dr (lower-right panel) for an overmetallic model ([Fe/H]=0.30) of 1.15 M and 1.782 Gyr, without over- 
shooting. For the echelle diagram, solid lines are for I — 0, dotted lines for I = 1, dashed lines for t = 2, and dotted-dashed lines 
for 1 = 3. 



where v is the frequency of the considered mode, r the radius, 
and c(r) the sound speeed. So that, for the 1-2 waves, the 
critical frequency for which the turning point is at the core edge 
is approximately 

v c ^0.4^. (5) 

Kcc 

Although c(R cc ) may slightly vary from model to model, it 
is not surprising to find that fixed values of v c corresponding to 
models with similar values of R cc . As a consequence, this be- 
haviour appears earlier for more massive stars and for models 
with similar sizes, but different luminosities, effective temper- 
atures, and gravities. 

3.4. Influence of the chemical composition. 

Tables [T] and |4] present the parameters of the transition models 
computed with the solar composition (here we use the Grevesse 
& Noels 1993 composition, compatible with helioseismology). 
As will be discussed below, Table [4] includes overshooting, 
which is not the case for Table Q] Tables [2] and [5] correspond 
to models computed with a larger abundance of metals ([Fe/H] 
= 0.30) and a helium value following the metallicity-helium re- 
lation obtained for the chemical evolution of galaxies (Izotov & 



Thuan 2004). Finally Tables|3]and|6]correspond to models with 
the same overmetallicity but with a solar helium abundance Y Q . 

The strongest effect of overmetallicity is to move the evo- 
lutionary tracks to lower effective temperatures. This is true for 
both helium values. In the case of the large helium mass frac- 
tion, the ages of the overmetallic transition models are sim- 
ilar to the case of solar values; however, this corresponds to 
different stages of evolution. For example, the models of 1.15 
M presented in Tables[TJand|2]have similar ages, but the one 
with overmetallicity and large helium abundance (Table 13 is 
still on the main sequence with a convective core, whereas the 
one with solar abundances (Table Q]) is already on the subgiant 
branch, with a helium core but no convective core anymore. As 
can be checked in the tables, their luminosities, effective tem- 
peratures, and log g are consequently quite different. 

For the case of overmetallic models with a solar helium 
abundance, the ages of the transition models are higher than 
for models with solar metallicity. As for the solar case, the less 
massive models are already on the subgiant branch with no con- 
vective core left, but their luminosities and effective tempera- 
tures are lower. 
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Fig. 6. Small separations (upper-left panel), echelle diagram (upper-right panel), sound speed profile (lower-left panel), and 
integral I(r) = £ \j- r dr (lower-right panel) for an overmetallic model ([Fe/H]=0.30) of 1.15 M and 4.654 Gyr, without over- 
shooting. It is the first model for which we obtain negative small separations below 3.5 mHz. This model has a convective core. 
There is an important discontinuity in its sound speed profile at the boundary of the core. For the echelle diagram, solid lines are 
for I = 0, dotted lines for I = 1, dashed lines for C = 2, and dotted-dashed lines for £ = 3. 



3.5. Influence of overshooting 

As discussed in Sect. 2, overshooting is added as an extension 
of the convective core, with a thickness equal to 0.20 times 
the pressure height scale. Figures [8] [9] and [TOl present the in- 
ternal parameters of stars along an evolutionary track of 1.15 
M , [Fe/H]=0.30 with overshooting. These figures are directly 
comparable to Figures [5] |6] and [7] which presented the same 
parameters for stars without overshooting. When overshooting 
is added, the radius of the mixed central core (convection plus 
overshoot) is increased so that the negative small separations 
appear for lower frequencies. 

The overshooting cases are presented in Tables[4]to[6] The 
transition models appear at similar ages as for the cases with- 
out overshooting. However, as can be seen in Fig. [2] this cor- 
responds to a relatively earlier stage on the main sequence, as 
the evolution time scales are longer in the case of larger mixed 
zones in the stellar centre. 

In the overshooting case, due to the larger core radius, 
the oscillations predicted for the small separations (Bazot & 
Vauclair 2004 and references therein) begin to appear clearly 
(Figs. l9land[T0b. For example, in the case of 1.15 M Q , 6.793 
Gyr (Fig. [Toll, the small separations for t — - 1 — 2 begin to 
decrease for increasing frequencies, become negative at v = 1.8 



mHz, go on decreasing, go through a minimum at 2.7 mHz, and 
then begin to increase. They eventually become positive again, 
for a frequency higher than 3.5 mHz. 

This behaviour is related to partial reflexion of the waves 
inside the mixed core, and the modulation period (in the fre- 
quency space) is of the order of the inverse of the acoustic time 
corresponding to the core boundary (given in the last column of 
Tables [TJ to [5). That we only observe this behaviour for models 
with overshooting is due to the increase in the core size leading 
to a decrease in the modulation period. 

4. Summary and discussion 

As powerful as asteroseismology may be, it will always give 
much less precise results about stellar internal structure than 
helioseismology does for the Sun. While more than ten mil- 
lion modes can be identified for the Sun, only tens of modes 
may be observed for stars. In this framework, skilled recipes 
are needed to better characterise the stellar structures from the 
observed acoustic frequencies. Among the frequency combi- 
nations that have been proposed to test stellar interiors, the 
"small separations", which represent the frequency differences 
between modes I, n, and i+2, n— 1, are known to give important 
hints about the stellar cores. 
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Fig. 7. Small separations (upper-left panel), echelle diagram (upper-right panel), sound speed profile (lower-left panel), and inte- 
gral I(r) — J r i j^dr (lower-right panel) for an overmetallic model ([Fe/H]=0.30) of 1 . 15 M and 5.937 Gyr without overshooting. 
This model is on the subgiant branch with a helium-rich core. Here the integral I(r) changes sign near the core. For the echelle 
diagram, solid lines are for i — 0, dotted lines for £ = 1, dashed lines for 1=2, and dotted-dashed lines for I = 3. 



In the present paper, we have shown that the small separa- 
tions, which should, in a first approximation, have a positive 
sign, may become negative during stellar evolution. Moreover, 
this behaviour occurs systematically for solar-type stars and 
may become observable at the end of the main-sequence or at 
the beginning of the subgiant branch for the t = - 1 = 2 sepa- 
rations. It also occurs for the I = 1 - € = 3 separations, at higher 
frequencies that are generally outside the observing range. For 
this reason, we concentrated on the t = - i = 2 ones. We 
found that this effect is directly related to the presence of a 
helium core, itself induced by convection in the central stel- 
lar regions. Although convective cores by themselves lead to a 
significant decrease in the sound velocities in the central stel- 
lar regions, the resulting gradients are not sufficient to inverse 
the sign of the small separations, until enough helium has accu- 
mulated. At the end of the main sequence, or at the beginning 
of the subgiant branch, the chemical gradient induces a much 
larger sound velocity gradient than convection alone and leads 
to negative small separations in the observing frequency region. 
When the star lies at the beginning of the subgiant branch, the 
convective core has disappeared but the effect continues due to 
the remaining helium core. From then on, this behaviour goes 
on, until the stellar structure becomes strongly modified on the 
subgiant branch. 



The occurrence of negative small separations stems from 
the modes of different I not seeing exactly the same stellar re- 
gions. In particular, while the I = modes travel down to the 
stellar centre, the I = 2 modes travel between the stellar surface 
and a "turning point" that depends on its frequency. Modes with 
low frequencies may penetrate the stellar core, whereas modes 
with high frequencies stay outside. The £ = 2 modes for which 
the turning point corresponds to the edge of the stellar core 
play a special role in this respect. In most cases, negative small 
separations begin to occur at the frequency of these modes. 

The interest of this behaviour of the small separations is 
evident, as it may help for characterising convective and he- 
lium cores and for constraining overshooting, at least for stars 
at the end of the main sequence or the beginning of the sub- 
giant branch. However, this may not be very simple in real life. 
When we observe the acoustic spectra of a solar-type star, the 
identification of the modes is not easy. That there may be a 
mode crossing point adds a new difficulty in this respect. Once 
an observed "echelle diagram" has been obtained, it may be 
difficult to decide whether it shows a crossing point or not. 

We will discuss these observational versus theoretical 
points in a forthcoming paper, in which we will give a new anal- 
ysis of the star // Arae (after Bazot et al. 12005). We will show 
that models with crossing points for the t = - I = 2 modes 
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Fig. 8. Small separations (upper-left panel), echelle diagram (upper-right panel), sound speed profile (lower-left panel), and inte- 
gral I(r) = £ \j^dr (lower-right panel) for an overmetallic model ([Fe/H]=0.30) of 1.15 M Q and 2.216 Gyr with overshooting. 
For the echelle diagram, solid lines are for i — 0, dotted lines for ( = 1, dashed lines for £ = 2, and dotted-dashed lines for i — 3. 



can be obtained if overshooting is added to the convective core. 
However, only very specific models may fit the overall echelle 
diagram, so that the constraints are really strong. We will also 
show that, in the presence of rotational splitting, the modes can 
be separated without ambiguity, as I — modes cannot show 
the splitting effects. 

As a conclusion, the negative small separations are quite a 
helpful tool for characterising convective and helium cores, and 
they may give strong constraints for the possible overshooting 
at the edge of the convective cores. 
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Fig. 9. Small separations (upper-left panel), echelle diagram (upper-right panel), sound speed profile (lower-left panel), and inte- 
gral /(r) = J r R \j- r dr (lower-right panel) for an overmetallic model ([Fe/H]=0.30) of 1.15 M G and 5.417 Gyr, with overshooting. 
It is the first model for which we obtain negative small separations below 3.5 mHz. For the echelle diagram, solid lines are for 
I = 0, dotted lines for I = 1, dashed lines for 1-2, and dotted-dashed lines for I = 3. 
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Fig. 10. Small separations (upper-left panel), echelle diagram (upper-right panel), sound speed profile (lower-left panel), and 
integral I(r) = ^ \%d r (lower-right panel) for an overmetallic model ([Fe/H]=0.30) of 1.15 M and 6.793 Gyr, without over- 
shooting. This model is on the subgiant branch, with a helium-rich core. For the echelle diagram, solid lines are for I — 0, dotted 
lines for ( = 1, dashed lines for ( = 2, and dotted-dashed lines for ( = 3. 



